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SUMMARY 



There is e3q]erlraental evidence that channel flows involving 
shock-free deceleration through the speed of sound are unstable. 

An analysis of ncnviscoue imsteady channel flow has- "been made to 
gain scans insist into this apparent instability, to provide some 
information on the factors determining the minimum shock intensity 
for a stable flow, and to study in general the formation afld motion 
of shocks in channel flo't'^ . It Is hoped that this treatment will 
be useful in two ways. First, the theory will have direct appli- 
cation to the study of normal shocks occiuring in supersonic nozzles 
and diff\xsers. Second, it might bo useful in leading to the 
solution of the more difficult problems of the formation and 
stability of shocks in two and throe dimensions. 

The analysis is divided into two parts. In pert I, Riemann's 
theory of tho propagation of finite-amplitude disturbances in a 
homogeneous medium is extended to tho case where upstream-moving 
pulses are superposed on a decelerating transonic channel flow. 

Shock waves are formed as the pulse is propagated, as in Riemann's 
problem. If the intensity of ttie shocks is essumed to be very small, 
the pulse approaches a ti-apped" state in which a portion of the , 
decelerating channel flow is converted to an accelerating flow which 
is an alternate steady-state solution for the channel. To the order 
of accuracy of the theory of part I, the trapped pulse becomes 
stationary in the sonic region of the channel and, thus, the flow 
has neutral stability. 

In part II, the assumption of weak shocks is dropped and the 
motion of sh'ock waves is considered more accurately. It is found 
that trapped expansion pulses are consumed by the shock motion and 
trapped compression pulses Inevitably grow. Thus, it is concluded 
that smooth transonic deceleration is unstable to compression pulses 
coming from the rear of the channel. 
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A similar, analysis is. mads -for the case vbere’ the shock wave 
forms a part oif tho equilihrlum channel flow. The shock position 
is found to he stable in diverging' Channels and unstable in 
converging channels. This analysis leads a 3 so to suggestions for 
the design and testing of supersonic diffusers. 


TOTRODUCTITO 


Experimental investigation of transonic channel flows has 
established a striking difference between the processes of 
acceleration and decei.eratlon through the speed of sovind. In a 
de Laval nozzle, a gas is accelerated to the velocity of sound in 
the converging (that is, converging downstream) part of the nozzle 
and is accelerated further to supersonic velocities in the diverging 
part. If the nozzle , 1 s shaped to avoid the condensation of 
coiipression waves, the acceleration throvigh tho speed of sound can 
be made very smoothly (without shock waves) . 

As far as is known, the reverse flow, smooth deceleration 
through the speed of sound in a channel flow, has never been 
observed. The experimental situation can be summed up briefly as 
follows; In tho process of starting the supersonic flow, a normal 
shock always forma ahead of tho converging part of the channel 
(supersonic diffusc-r) . By suitably changing the velocities and 
pressures at the ends of the diffuser, the noi'mal shock can bo made 
to ilump to the diverging part. With further changes in velocities, 
pressiures, and/or geometry, the normal shock can be pushed up tho 
diverging channel and its intensity reduced some'^diat. However, 
before the shock can be made to disappear, it suddenly (to visual 
observation) Jumps to a position ahead of the converging part of the 
diffuser and the starting process must begin again. The experiments 
strongly indicate, therefore, -that smooth deceleration throu^ the 
speed of sound in a channel flow is unstable. 

An analysis of this apparent ins^t/ability is tho central 
problem of this paper. Specifically, the stability of transonic 
decelerating channel flow to a special type of disturbance coming 
frem the rear of the channel is considered. The analysis is 
restricted to channels of slowly varying cross section so ■Jhat a 
quasi cjne-dimenslonal treaianent can be used. This treatment is 
based on Rleraann^s treatment (reference l) of plane waves superposed 
on a homogeneous raodlum. A part of Eiemann's paper will be 
paraphrased herein to put his results into a form suitable for 
extension, to channel-flow problems. 

The author is very grateftil for the helpful discussion of 
these problems with Mr. A. Kahane and Mr. Lester Lees. 
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EIEMMn 'S TREA1MSCT OF THE PKOPAGATICN OF FINITE PLANE 
DISTURBANCES IN A HOMOGH^EOUS i^EDIUM 


In 1859 , Riemann puLliehed his great paper (reference l) on 
the propagation of plane aerigtl. wavea of finite amplitude. Starting 
vith the one -dimensional equations of motion and continuity and 
ass{mlng Isentropic flow, he shows (among other things) that: 

( 1 ) a disturhance initially confined to a finite region will spread 
out into waves propagated in hoth directions froa the soui’ce, and 

(2) that coExpreeslon fronts in the propagated waves will steepen 
Indefinitely, and will eventually form a ccmpression shock. 

Eietnam starts with the one -dimensional equations of motion 
and continuity (assuming isentropic flow of a perfect gas) : 


+ u = 0 ( 1 ) 

ot ex bx 


and 


* 0 (2) 

ot ax 


(Symbols are defined in appendix A.) 

He then multiplies equation (2) by and adds to equation (l) 

obtaining 


^ + a + (u + a) f ^ ± a * 0 (3) 

ot ot \,3x 8x / 

Assuming all changes to be isentropic and introducing 7 = l.^^ 
(the value for room- temperature air) yields 

d log p <= .. - s (h) 

(7 - l)a a 

Hence, a d log p = d 5a* Using this resxilt permits equation ( 3 ) 
to be written 


^ (u i 5a) + (u ± a) ^ (u t 5a) <= 0 


¥ 


(5) 



h 


HACA TN No, 1225 


It is convenient to introduce the syrnhols P = u + 5a 
and Q a u - 5a which occur frequently in these equations . 

With these syTibola u + a = and u - a « 2P_+__22^ 

5 5 

The equations (5) can he written 

and 


M + gP + ^Q. 

3t 5 hx 


0 


(7) 


If the dlsturhances are not too large (| u 1 < a) these 
equations represent dlsturhances moving in the ±x directions, 



If, for example, there are no dlstAirhances in a given region 
moving in the -x directicn, then Q ■= Const. = “5Qg where 

is the velocity of sound in air at rest in this region. In this 
case, it can he seen from equation (7) that Q remains constant, 

since ^ “ 0, even in the presence -of a P wave. Thus, waves 

moving in one direction in a homogeneous medium do nob reflect 
waves moving in the opposite direction. 


In the analysis of channel flow, finite pulses coming from 
the rear of the channel will he assumed. The term ’finite pulse’ 
is used to denote a finite disturbed region hounded on both sides 
by undisturbed fluid (or equilihrlua flow in a channel-flow 
problem) . A finite pulse will consist of at least one phase whore 
the. pressure drops, which vrlll be called the expansion phase, and 
at least caae ccmprepsion phase where the pressure rises with time. 

P P 

It can be shown that the integral ^ P dx (or J Q dx)’, which 

will be called the pulse area, taken over the whole of a finite 
pulse does not change with time ■vdlen pulses moving in only one 
direction are present. Integrating equation (6) with respect to i 
gives 

^y’p dx +J*' dP « 0 
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Since P vanishes at "both en(te of the pulse, the second integral 
is 0 for Q = Const, and j "B 6 lX (or similarly dx) is 
constant in time^. 

Eiemann used equations (6) and (7) to show that compression 
fronts in these waves steepen to form shocks while expansion waves 
flatten out indefinitely. This Infomation can he obtained for 
the case where only waves moving in one direction are present hy 
a different approach which will he readily extended to the problems 
considered later hy examining equation (7) • Differentiating 
equation (7) with respect to x, assuming P = Const., and intro- 
ducing e = there is obtained 

ox 


|£+ 2P+-2fi - 2,^ 

3t 5 5i 5 


( 8 ) 


For the case where only Q disturbances a3?e present the quantity € 
may be understood physically as twice the velocity gradient in the 

disttjrbance since u =« 2~L-£. 

2 


Consider now that the point under observation moves with the 

propagation velocity and let ^ denoto the rato of change of e 

dt 

for this point. Now 


§1 = 8® + 8x 
dt 5t dt §x 


where §5 is the local propagation velocity 3Q, Hence, 

cL'b ^ 

more general form of this . theorem can be obtained as follows; 
let f(p) be an arbitrary well-behaved function of P. Multiply 

equation (6) by integrate as before. It follows 

that f(P)dx = 0 where the integration is over the whole of 

a finite ptilse . Thus, the mass, impulse or energy carried by a 
pulse will also remain constant and any of these qxxantities could 
be used. The quantity ^/’Pdx Is convenient since P appears 

directly in the equations. 
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equation ( 8 ) can be -written 

at 5 

or 


£ = 1 

e ate + 1 

5 


where S is the value of c at t - 0 . Now, If £ Is positive, e 
decreases toward zero as -the wave Is propaga-ted. If S Is 

negative, |€( Increases Indefinitely, hecoirilng Infinite at t s= - 

3 ^ 

For the case where only Q waves are present In the region under 
consideration, positive and negative values of e can he Identified 
from -the fundamental equations as expansion and compression waves, 
respectively. Thus, compression waves steepen to form shocks and 
rarefaction waves flatten out indefinitely. 


GENERAL PLAN OF THE CHANNEL-FLOW ANALYSIS 


The general plan of this -ferea-tment can host he explained hy 
considering -the central problem, the stability of smooth deceleration 
•through the speod of sound, from a physical point of view. In such 
a flow, disturbances starting In -tiie rear of the channel will be 
trapped in the sonic region. The presence of -this disturbance trap 
constitutes a distinct difference between this flow and smoo-th flows 
kno^^n to be stable . 

Two typos of processes -vdiich govern -the behavior of pulses in 
the sonic region are considered separately in parts I and II of 
this papei*. First, tlaere are ism-tropic processes similar to the 
propagation of waves of finite amplitude considered by Rierannn. 

Shocks form in the, channel flow as in Eiemann's case. If -these 
shock waves are assumed to be of very small in -tensity, it is sho-wn 
in part I -that a pulse approaches a stationary s-tato. Thus, the 
amplitude, shape, end ;posltion of a pulse approach assyroptotlcally 
what will bo called a ^’trapped pulse.” The motion of the shock 
waves in these trapped pulses is considered more accurately in 
part II, and it is found that -the shock moves and causes -the piilse 
to be consumed or to grow, depending on its sign. In. cases where 
the intensity of the shocks involved is not -too large, the 
consumption process is much slower than -the trapping process . In 
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these cases, it is possible to obtain a fairly complete history of 
the behavior of a pulse without considering the interaction of 
these two processes. Thus, in part I, the rapid trapping process 
is considered and the slow consumption of the pulse by shocks is 
neglected. In part II, the consumption of fully trapped pulses 
by shock waves is discussed. 

A similar analysis is made for the case where a shock wave 
forms part of the initi.al channel flow. In this case, disturbances 
from the rear serve to move the shock to a new position and, if 
changes in shock sti'ength are neglected, it is found that the 
pulses in this case are also trapped as shock displacements (shown 
in part I). Here again, the slower effects on this shook displace- 
ment are discussed in part II, the changes in shock strengtli being 
included. In this way, the stability of shock waves in channel 
flows is studied. 

The boundary— layer effects are neglected entirely in this 
paper. It 1s well known, hov/ever, that such effects are important 
in channel flows involving shock waves. . The results obtained can 
therefore, at best, be only qualitatively similar to experimental 
results. Perhaps it will be possible eventually to add boundary- 
layer effects for cases vhore they a’-e Important, thus following 
a procedure frequently used in gas dynamics. 

It should be pointed out that ^nf ormti cn of the type sought 
in this analysis cannot be foimd by introducing the acoustic 
approxlmationj that is, that the disturbance velocities are negligible 
compared with the propagation velocity. In the sonic region of a 
channel, the propagation velocity of an upstream-moving pulse is 
close to zero. On the other hand, if the acoustic approximation is 
introduced, the disturbance velocities tend to become infinite in 
the sonic region (reference 2) . Thus, it is clear that the acoustic 
approximation cannot be used to study the development of upstream- 
moving pulses in the sonjc region. 


I - TEE TRAPPING OF FOLSES AND THE FORMATION OF SHOCKS 
IN A DECELERATING TRANSONIC CHAIiNEL FLOW 
The Propagation of Finite Plane Disturbances 
Superposed on a Channel Flow 


General equation s.— An analysis similar to Riemann’s can be 
applied to the case where plane disturbances are superposed on a 
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one-dliriensional flow field. (A oue-dimensional flow field is one 
where the variations in velocity, density, etc., perpendicular to 
the flow direction, are nesligihle compared with those in the flow 
direction.) Consider the field of flow in a straight channel of 
slowly varying cross section S. The eqmtion of continuity hecomes 


S^P ^ 8(PitS) 
8t 


= 0 


or 

^ Ipg- P + + ^ + = 0 (9) 

dt ox 8x dx 

Multinlying equation (9) hy ta and adding to equation (l) yields 


A(u ± 5a) + (u ± a)^(u ± 5a) ± ua^-iS^ = 0 (lO) 
dt ox dx 


An equilibrium (steady) flow in the channel with quantities denoted 
by subscrint o will be aosuned. Priiaod quantities will denote 
departures frora this equilibrium. Thus, 


( 11 ) 


u = Uq + u' *1^ 

a » + a' J 

By substitution of equations (ll) in equation (lO) , there is obtained 


^ 'u' t 5a') + (uo ± ao)”(uo t 5ao) * Upao ^ - ^ - ^- 




dx 


+ I (u„ + ap) + tu’ ± a')l L^a J I + (u> + an<l.( uo. - -^ap) 

o Q • J dx dx 


i (u'a' 


+ u'a + a 'u)- *^ ® 0 

o ' dx 


( 12 ) 
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Note that, where ^ i 


a02g^ . - 1)!^ 

dx u_ ^ dx 


o 


and 


dx ^ dx 


Dropping the sec end and third terns of eqiiatlc®. (l2 ) , wfaich cancel 
each other (equation (lO) for steady-flow case), gives 


|j{u' t 5a’) * [fu„ ± ao) (»' ±.a')J 

. fu' i a')U jMo)^ i (Wo® - * ^.+ . 0 ( 15 ) 

dx ^ \ % Mjj / dx 

If the quantities P and Q are introduced, where now 
P = u' + 5a' and Q « u' - 5a’, the equations ( 13 ) becoiaei 


db V' 


+ a,^ 
o o 




P - c 


^•10 

) dx 




and 





P + Q. + P - qN^^^o 
2M^ 10 / dx 


0 


(15) 
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Reflections In a velocity field .- The equations (l^^■) and (15) 
indicate that disturbances will ho propagated upstream and down- 
stream in the channel flow as was the case in Riemann’s problem. 
Because of the last two terms in these equations, however, a pulse 
moving in one direction will reflect disturbances moving in the 
other direction. Perhaps it would be possible to solve these 
equations numerically for any particular case by the method of 
characteristics. It will, however, serve the general purposes of 
this paper better to derive general approximations rather than more 
accurate numerical results for spoclfic cases. 

In part I of this treatment, the reflections will be neglected; 
but before this assumption is introduced, an evaluation of the order 
of magnitude of these reflections will be made to provide an 
indication of the error thus introduced. 


First, consider a pulse moving upstream into an equilibrlClm 
velocity field; that is, containing no downstreamr-movlng disturbances. 
Consider a point moving with the downstream velocity of propagation 

u + a * Uq + ao + — -E— The rate of change of P at this 

P SP 

observation point will be denoted by ~ and is equal to 

^ rate of growth of P while this point crosses the Q 

ot dt ox 

pulse is thus given by the STim of the first two terms of equation (l4) 
that is 



(1 - Mo) 


3P + 2Q 



^ P + Q 

2Mo 



In order to obtain a simple estimate of the reflections produced in 
the sonic region of a channel, it will be convenient to introduce 
the following approximations in the right-hand side of this equation; 
(1) P is negligible compared with Q, (2) Mq + 1 = 2Mq = 2, 

, . q2 q 

(3) 25u^ negligible compared with g. Thus, there is obtained 


dt 



Q |(1 - M^) 
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Now, integrating over the >41016 of an upstream moving pulse gives 
for the value of P at the rear (which will he the maximum value 
of P if Q is all of one sign througjiout the pulse) 


^rear “ 5 X 

^pulse 

Since thfe pronagation velocity is close to Sag, dt 
■Doint under ohaervation) , and approxi-mtely 


(for tte 

2^0 


P 

rear 


2 

5 


(1 


%uil 


se 


- 

2^0 



dx 


For example, for a square pulse - that is, a pulse for which Q = Const, 
from a point where the steady velocity is u^ to a point where it 
is Uq + AUq, there is obtained 


P = Q 

rear 




^o ^ 


i-d - Mo) 


av 


where (l - I'io)av average over the region of the pulse. 


The relative magnitudes of ? and Q in this problem can be 
Illustrated by a n’mnerical example. If (l - = 0.1 


Auo 

and 0 .1, 

^o 


then P 


rear 


-= - -SL, 
500 


It thus apnears that for short upstream-moving pulses in the 
sonic region, it will be possible to neglect P in comparison 
with Q in equation (15) if there are no downs treem-rioving 
disturbances ureBent Initially. 

Conservation of -pulse area .- If reflections are neglected, the 
propagation of an unstrea:i-moving- or Q pulse can be studied with 
the aid of equation (15) alone. Since this treatiaent is primarily 
concerned with the sonic region of charnels, - it wi3JL be possible 
to obtain si:mle apmroxi-riations by neglecting the fourth term of 
equation ( 15 )^ >7-hich contains the factor .(Mq^ “ l)^ comparison 
with the third term. Eliminating this term, replacing 1 + Mq by 2 
in the third terra, and neglecting P, compared with Q, gives 
for equation (I 5 ) 
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if*(uo-=c*f|S.|,?g = 0 ( 16 ) 

From Riemarin's equations, it was shown that the quantity J Q dx 

(pulse area) is constant in the isentropic propagation of pulses 
in a homogeneous gas . By the use of the approximate equation (l6) , 
it can easily he shovn that the same rule Is approximately true for 
propagation in a flow field. Equation (l6) can also he written 


St 55 


I * “o - ^o) 


= 0 


(17) 


as can readily he verified hy differentiation (for values of a 
close to l) . 


Integrating this equation with respect to x over the range x^ 
to x-^ gives 



If the integration Is taken over the whole of a finite pulse, the 
second term vanishes and the pulse area J Qdx remains constant to 

the order of accuracy of equation (l6) . It can also he seen hy 
taking the point x^^ ahead of the pulse where Q » 0 in equation (l8) 

that, at an arhltrary point x^, Q pulse area flows through at 
the rate f-Q/^ fi + u_ - a^V"} 

L 2 ° V.U 

It can he shovn that, as long as on3.y upstream-moving disturhances 
are' present, this rule - conservation of pulse area in upstream- 
moving pulses - is not affected hy the presence of shocks (if do^stream 
moving disturhances produced hy the shocks are neglected) . For the 
case where the shock is not part of the initial equillhrlum flow (that 
is where it develops in an upstream-moving distur’bance) and P « 0, 
it is shown in appendix B (equation (BV)) that the velocity U of 
propagation of a shock wave is 


U = Uq - a^ + ^ 


(19) 
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•where -raluea of Q before after the shock are denoted by 
0^ and Qg, respectively. 


If now it can be shown that the pifLse area entering the shock 
per unit tine is equal to that lea-ving the shock, the pulse 
conser-viation rule will be satisfied. The pulse area entering or 
leaving the shock is the sum of the pulse area passing throu^ a 
B-tationary point and the pulse area overtaken by the shock j thus, 
it must be shown that 

UQi - ^ + Uq - ao) = UQg - ^ + Uq - (20) 

If the expression (19) le substi'bu'ted for U, this equation 
is identically satisfied, and thiis, this extension of -the rule is 
established. 

Another similar case which will be of interest in connection 
with later considerations is that for an upstream-moving pulse of 
small amplitude interacting wj.-th a shock that is part of the initial 
equilibrium flow. This case is illtistratod in figure 1. It is 
necessary to assume that the shock intensity is low (equilibrium-shock 
Mach number close to 1) bo get a simple result. It is also necessary 
to neglect the variations wi-th x of the local equilibria velocities 
and velocities of soundj in comparison with the disturbance velocities. 
These changes, which eire due to the changing channel area, are 
responsible for the stability of position of the shock wave. Thus,, 
in other words, the effects due to the stability of the shock wave 
are neglected as compared ■with the dist'urbance effects. The 
circumstances under which this assui^ption is permissible will become 
cleeirer frcsm the results of part II. If the stability of a normal 
shock is neglected, the velocity with which it moves tinder -the 
influence of a small upstream-moving disturbance is calculated in 
appendix B (equation (B9)) to be simply 

U = 

In addition to the contributions to the pulse area which were 
present in the previous problem, pulse area is con'tributed by the 
fact that the shock is merely out of position. The equation which 
must now be satisfied to demonstrate conservation of pulse area is 
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vhere tiie subscripts find 2^ denote equilibrium conditions 

before and after the shock, respectively. Noting that approxi- 
mately (see appendix B) 


and 


Uo - ag , 

Ui - ug “ 5 

•^o o 



for Mach numbers close to 1 and introducing the shock velocity 
(equation (21)), satisfies equation (22) identically. Thus, the 
rule of conservation of pulse area anpears to be a usefixl approxi- 
mation with considerable generality. 

The nulse shaue The shape of a pulse (that is, Q as a 
function of x) advancing upstream in a velocity field, in the 
absence of dovnstrean-moving disturbances, can be studied in a flow 

field by considering its x derivative s g . The variation of 

this quantity with time can be found by differentiating equation (l6)‘ 
with respect to x. 




5 \ 



For the case where the equilibrium velocity gradient is constant, 
cL^y. 

that is — -2 = 0, the behavior of the pulse shape can be studied 
dxr 

simply from this differential equation. Consider the point under 
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oliservation to move with the wave, that is, move with the local 

propagation velocity. Let ^ he the rate of change of e when 

dt 

the noint mder ohservatlon moves with the propagation velocity} 
■i±en the preceding equation hecomes 


M. 

dt 


St \ c 


+ 1 
5 



(23) 


The • variables in equation ( 23 ) can he separated. If the 


dimensionless parameters E 


€ 


and T 


4 . L- 2 . 

dx du 

equation ( 23 ) can he integrated for — 2. 

dx 




dx 

constant to give 


are introduced, 


E 





where f is the dhriensionless slope of the wave at the time t = 0. 
Positive values of T in equatim (24-) will represent accelerating 

equilihrium-velocity fields y — 2> 0 | and negative values of T 

will represent decelerating fields. Also, in accei.erating fields, 
positive values of E will represent expansions and negative 
values will represent compressions, end the opposite is true for 
decolcrating f ialds^ Equation (24-) is plotted in figure 2 for 
various values of E. 

It can he seen from this figure that the results are considerahly 
different from Riemann’s. Thus, in an accelerating field, compression 
waves must have a certain initial steepness, S < -1, or they will 
flatten out in time. On the other hand, in a decelerating field, 
expansion waves, instead of flattening out (as occurs in Riemann’s 
case), apnroach a slope given hy S *s -l. It should he noted 

that E ts -1 coj.’resnonda to a tota3. velocity gradient ^ - £122.. 

dx di 

"Tranued" nulsee .- The approximat® relations which have been 
developed In the preceding sections can he used to show that a 
nulse moving upstream in a uniformly decelerating flow approaches 
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a Btationary state. Thu^, the pulse area remains constant, the 
expansion phases approach the slope given hy E = -1, and the 
compression phases ‘become shocks . As these developmen+'S take 
place, the piUse advances to a position wheivt) the propagation 
velocities of the shocks approach ^ero. tflien^^this position is 
reached the pulse vill be said to be trapped* • It should be pointed 
out that the trapped pu3.se is an idealization and would not exist 
in an exact treatment. The introduction of this concept makes 
it possible to divide the history of a nulse into two periods 
which can be considered separately. The trapping of compression 
and expansion pulses in a decelerating transonic channel flow is 
shown schema tlcal3.y in figwe 3* 

The asymptotic result, for expansion phases of trapped pulses 
E = -1 is co.mpared in figure 4 with the exact accelerating 
equli.lbrlum channel flow. It can h© seen from this figure that if 
all the velocities Involved are close to the local speed of sound, 
the expansion phases are very close to the accelerating steady flow. 

It seems reasonable to infer from, figure 4 that the differences 
between the two curves are due to the approximations introduced in 
deriving equation (l6) from the exact equations (l4) and (15) 

(neglect of reflections and - 1 terras) . In other words, it 

seems reasonable to infer that an exact calculation from equations (l4) 
and ( 15 ) would yield the exact accelerating steady flow. This 
inference will be adopted in part II and the exact accelerating flew 
will be used in place of the E = -1 result in caD-culations of the 
shock velocities . 


It will be valuable to calciilate the order of magnitude of the 
times in which the first-order trapping processes take place. These 
times will be compared with times in which the slower procBsses of 
part II take place. For example, the leading edge of an expansion 


wave propagates upstream with the velocity u. 


*■ 0 * 


If X is the 


distance of -tois leading edge from the sonic point of a channel, 
then In a miformly decelerating field (for close to 1) 

6 

u„ - a„ a -urtr x . The equation for the approach of the leading 


o.x 


edge to I 'he sonic point is thus. 


dx 6. 
dt 5 dx 


The relaxation time for the approach, is therefore 
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A relaxation, time for the slope of an expansion wave to approach 
its equilibrium slope can similarly be found from, equation (23) • 

d.u^ 

If the slope is close to the equilibrium value e = ~dx^ 
equation ( 23 ) can be written 



+ 4 



) 


The relaxation time for this process is then or one half 

5 V to/ 

as large as that for the approach of the leading edge to the sonic 
point. It is clear frm figure 2 that the formation of shocks from 
compression waves is a faster process and a shock is formed in a 
finite time. 


Application to the Formation of a Normal Shock 
in a de Laval Nozzle 

The theory Just developed can be applied to trace the formation 
of shock waves in a de Laval nozzle. Consider the converging- 
diverging nozzle shown schematically in figure 5* Suppose the back 
nresGure to be adjusted to such a value that the local speed of 
sound is reached at the throat but not exceeded. The form of 
velocity distribution that would be obtained is also shown schemati- 
cally in figure 5* The small curved region near the local speed of 
sound in the velocity curve which -would be obtained in a practical 
case has been neglected to sim-olify the argument. 

It is kno-wn frcxa mejiy experimental results -that if the back 
pressure is levered further, a normal shock will form in the nozzle. 
The fo3riation of this shock can be readily described by use of -the 
foregoing theory if the velocity at the end of the nozzle is 
sufficiently close to the speed of sound. It will be assumed -that 
the back pressure will be lowered continuously for a time and then 
held constant, (it will al.so be assumed that the back pressure 
is not lowered enovi^ -to produce supersonic flow at the end of the 
nozzle.) While the back pressure is being lowered, the velocity 
at the end of the nozzle will increase . The space derivative of 
the velocity disturbance created can be found frem equation (l6) 
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If downs tream-moving disturbances are neglected. Solving this 

equation for ^ ^ yields 

dx ? dx 


dx 


ft 


6 

1 


u 




u' 



du’ 





(25) 


It can be seen from equation (25) ■ttiat some time before the expansion 
stops, will become negative. This result will occur when 


6 

5 


u 


^ 0 
dx 5t 


since the denominators In equation (25) are negative for subsonic 
velocities . 


A negative value of means that a compression wave has been 

produced in the uroceas of decreasing the back pressure to a new 
steady value. (It shoxtLd be remarked that no such catipression waves 
are started in the absence of a decelerating velocity field.) This 
compression wave will travel upstream and will steepen as has been 
shown previously (equation (23)) to form a cnrapression shock. The 
compression shock and the wedge •'shaped, area ahead of it will grow 
by the addition of expansion pulse area coming from the rear of the 
nozzle (pulse area is continually contributed since Q / 0 at the 
rear) . The shook will grow until its intensity end its dissipation 
are large enough to reach a condition where the nozzle flow is 
again at equilibrium. 

The oalculations on which figure 5 is based are given in 
appendix C . The leading phase which originated while the nozzle 
back press\3re was falling can be calculated from equation (23) and 
•the laiov/n propagation velocity of its leading edge. The s’beady 
phase which originates after -the nozzle back pressure has again 
become steady can be found by ii^tegratlng equation (l6), use 

being made of the fact that 2JL. a o. for •this phase. For the 

ot 

case considered, an infinite decelerating velocity gradient appears 
at the point where the velocity first reaches the local speed of 
sound. The shock that follows is calculated s'bep by step frcsu 
equation (l9) • The calculations used in plotting figure 5 ar© thus 
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only a first approximation to the shock motion. A more accurate 
shock velocity calculation (by the methods of part II) would show 
the shock to settle at an equilibrium position. 

If the unsteady aspect of a channel-flow problem is considered, 
it is possible to see how shocks arise in the transonic flow. It 
is quite possible that an analysis of this kind could be extended 
to two or three dimensions and to the flow over obstacles and that 
if these extensions were made, the problem of the formation of 
normal shocks in these more complicated cases could be greatly 
clarified. 


II - THE MOTICN MD STABILITY OF CHOCK WAVES 
IN A CHi®HEL FLOW 
Preliminary Considerations 


The trapping of pulses and the formation of shock waves in 
decelerating channel flows have been considered in part I. In the 
case of the central problem, the stability of smooth transonic 
deceleration, these calculations indicate that a trapped pulse 
converts a nart of the channel flow to the accelerating steady flow. 

The trapped pulses always Involve shock waves and if the approxi- 
mations of appendix B are used for the shock velocity and if down- 
stream waves created by the shocks are neglected, it is found that 
the entire pulse becomes statjonary in the sonic region of the 
channel . 

The occurrence of these trapped pulses in transonic decelerating 
channel flow means that, to the order of accuracy of these calculations, 
this flow has a kind of neutral stability to pulses ccaning from the 
rear of the channel. Therefore, it will be necessary to make more 
accurate shock-velocity calculations before it can be decided whether 
this flow is stable or unstable. In this part of the paper, the 
motion of shock waves will be considered more accurately. 

It wi].l be assumed that conditions at the ends of the channel 
are steady; that is, the shock velocity will be calculated for the 
case where the channel flow is out of equilibrium (for example, vhere 
a chaimel shock has been displaced by a pul.se) but where no new 
disturbances are originating at the end of the channel. Thus, 
disturbances will be treated by the methods of part I and these 
more cumbersome calculations will be used to solve stability problems 
which cannot be solved by the methods of part I. 
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The flow on the upstream side of the shock will not "be related 
in any way to the shock position or motion hecause any disturbances 
created by the shock will not move, upstream as rapidly as the shock 
does. It is therefore possible to calculate the state of the air 
at all points ahead of the shock wave before the shock mption is 
determined. For all the applications of the shock- velocity calculations 
that are contemplatpd In this paper, an equilibrium flow (either the 
original equilibrium flow or the alternate one inferred to be produced 
by a trapped pulse) will be assumed upstream of the shock wave. 

In steady flow, once the conditions ahead of the shock are 
knot/n, it is possible to specify immediately from the Bankine-Hugoniot 
relatione the conditions behind the shock. In unsteady flow, however, 
since the shock velocity adds another unknown to the problem, it will 
be necessary to have another relation 3 n addition to the Eankine-Hugoniot 
relations before the shock Velocity and the conditions behind the 
shock can be specified. Since the conditions upstream from the shock 
wave are already fully determined, this additional Information will 
be a relation between the variables of condition downstream from the 
shock. The nature of this relation can be brought out most .simply 
by considering a particular case, namely, the case of a channel 
discharging into the ataosphere. It will be accurate enough to say 
that discharging into the constant-pressure atmosphere is equivalent 
to a constant pressure boundary condition at the end of the nozzle. 

If a P ■ disturbance moves downstream in the charnel'’, it will of 
course, involve pressure changes and, when the end of the nozzle is . 
reached, this disturbance will bp reflected as an upstream-moving 
disturbance. Thus, it will be possible for a given P disturbance 
to calculate the Q disturbance at the end of the nozzle. Now, 
applying the isentropic wave propagation equations permits the 
extension of this informatics to regions upstream. Thus, in 
general, if the downstream-moving disturbances at a point are known 
and the boundary conditions at the end of a channel are specified, 
it wouil.d be possible to calculate the upstream-moving disturbances 
at the same point in the channel. More accxu-ately stated, the 
downstream-moving disturbances for a previous- interval must be known 
in order to calculate the upstream-moving disturbances . 

These problems could probably be treated accurately by the 
method of characteristics and calculations of this kind would certainly 
be desirable. However, for the exploratory purposes of this paper, 
an approximation will be introduced to make lb possible to reduce 
these calculations to a closed form. Thus, the reflection conditicn 
will be appljed immediately behind the shock instead of at the end 
of the channel. For example, in treating the constant-pressure end 
condition, the pressure immediately behind the shock will, be assumed 
to be the subsonic equilibrium pressure appropriate to the shock 
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position. The reflection condition, in combination with the 
conservation laws, may then be used to obtain an eilgebraic expiresslon 
for the shock velocities. This procedure introduces the following 
errors; First, the changes in the strength of the waves in moving 
downstream to -the channel end and in moVing upstream to the shock 
are neglected. Second, the time required for a wave to move from 
the shock to the" channel end and for the reflection to move back 
to the shock has been neglected. In cases where the channel 
following the shock is not too long, this effect will not be serious 
for shock Mach numbers close to 1, since the shock velocities 
found are s»iiall compared with the local propagation velocities . 


Reflection Conditions 

The analysis will be made for three types of reflection 
conditions; constant velocity at the channel end, no reflections, 
and constant pressure. These cases will provide a considerable 
variety of reflections; thus a method is provided for estimating 
the slgaiflcance of reflections in these problems and, hence, the 
significance of the two types of errors Just mentioned. A further 
discussion of these errors will be given after the shock- velocity 
calculations have been presented. 

Constant- velocity c h annel exit .- Constant velocity at the 
channel exit could be approached practically by arranging a channel 
to discharge into a machine which takes in a constant volume per 
unit time. If, in such a setup, a disturbance moves down the channel 
from the shock wave, it will be reflected from, the end of the channel 
in much the same way as a sound wave is reflected from a closed 
organ-pipe end. The velocity will thus remain constant at the channel 
end and, in this case, u' = 0 at the end of the channel end P = -Q 
at that point. As discussed previously, this boundary condition will 
actually be applied immediately behind the shock. 

The shock relation is introduced now in the form 


^ 1^2 



e 




( 26 ) 


where and a^^ are the velocity and velocity of sound immediately 

upstream of the shock in coordinates fixed with respect to. the shock 
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wave and V2 and ag ere the velocity and velocity of sound 

Iriimedlately following the shock In these coorddnates . The notation 
to ho used in connection with problems of the motion of shock waves 
is illustrated in figure 6. Peturning to a coordinate system fired 
with respect to the channel yields for equation (26) 


- u)^ - u) = - vf * 

» / Ug - vj + 50 / (27) 


The constant-velocity reflection condition yields the information 
u'g SB 0 and hence Ug = Using this information and eliminating 

ag between the equation (27) gives the following solution for the 
shock velocity. 


U 


5 





“ 30Uj^ U£ + 25 ag® 

'i 'i 


m 


It shovild ])e poi:t^ted out that the plus sign corresponds to the 
expansion shock’ which violates the second law of thermodynamics 
and, therefore, the minus sign must be used. 


Wonreflectlng channel exit .- The condition where no waves are 
reflected from the rear of the channel could bo approached in the 
case where the region of the chsinnel containing the shock wave 
under consideration is followed by a sonic throat and a supersonic 
region. In this case, the only reflections present would be those 
produced by the variations in channel area between the shock and 
the throat. If all the Mach numbers are close to unity, the 
variation in channel area would be small and, therefore, the reflections 
should be email. 


The absence of reflections means tliat no upstream-moving 
distincbance will exist behind the shock in the absence of external 
disturbances. The reflection condition in this case is then Qg = 0. 

Applying this relation immediately behind the shock gives 

Qq a u*2 " 5a’2 ^^sing this Information in equations (27) makes 

it again possible to oli:!ilnate u'2 and a ’2 and to obtain 

+ AU-'^ + BI]^ + CU + D « 0 


(29) 
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vhere 
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Here again. It is possible to solve for IT If Uv eii<i a, are 

knovn . Of . the four roots of equation (29) , two are Imaginary and 
one corresponds to the expansion shock so that only one has 
physical significance. 

Constant-pressure channel exit .- The third boundary condition 
which will be considered is the case of constant-pressure channel 
exit. This condition could be apuroached in a channel that 
discharges into the atmosphere # Transferring this boundary condition 
from the channel end to the shock position gives the information 
that the uressure behind the shock is the same as the pressure which 
would exist at that point in the steady flow. The constant-pressure 
boundary condition can thus be stated as pg » Using the 

shock relation 


(% - u)^ ■ 1 

- 2 6 


(30) 


gives an immediate solution for the shock velocity . 
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The Stability of Smooth Deceleratlcai throu^ 
the Speed of Sound 

In the case of a chaimel flov involving smooth deceleration 
through the speed of sound, it was found in part I that a pulse 
coming from the rear of the channel was trapped in the sonic 
region of tlie channel. Tills result means that, to the order of 
accuracy of these calculations, the flow was neutrally stable to 
pulses coming from the rear and that more accurate calculations 
are needed to decide whether this flow is stable or unstable . 

Two kinds of approximations were made in deriving the trapped- 
pulse result. First, the reflections produced by isentroplc waves 
and the terms containing - 1 as a factor in equations (lU) 

and (15) were neglected. The expansion waves were found to replace 
a part of the original decelerating equilibrium flow with a flow 
close to the equilibrium accelerating flow which is also possible 
for this channel. (See fig. 4 .) It was then Inferred that exact 
calculations would yield the exact accelerating equilibrium flow. 

If this Inference la correct, the exact final state of an expansion 
wave has been obtained by this method so that the approximations 
introduced in the calculations of the escpansion waves will have no 
effect on the final results. 

Second, rough approximations to the shock velocity developed 
in appendix B were used and, with these approximations, it was 
fomd- that the shock neither consumed the pulce nor made it grow. 

Hence, any change in the shock velocity from these approximate 
values would lead to consumption or growth of the trapped pulses 
and, thus, the more accvirate shock-velocity formulas Just derived 
can be used to determine the stability of smooth deceleration 
through the speed of sound. 

For the case where the equilibrium flow is smooth deceleration 
through the epesd of sound, the quantities appearing in eqiiations ( 28 ) 
through (30) can I'eadlly be evaluated and the shock velocity calculated 
It will be seen from these equations that the shock velocity depends 
only on local Oondlbions at the shock position. These local 
conditions, of course, are all functions of a single parameter and 
the parameter used to present these results is the Mach number 
Immediately ahead of the shock. 

Results found for the three reflection conditions are plotted 
In figure 7 « It will be seen that the shock velocity is always 
negative; that is, the shock always moves upstream. Thus, in the 
case of the trapped expansion pulse (fig. 3) > the shock will move 
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UTDstream and eventually will eonsuaie the pulse. Smooth deceleratlcn. 
through the speed of, sound is therefore stable to exoansion pulses 
coming from the rear. On the other hand, in the case of the trauped 
compression pulse, the shock wave ■will again mo've upstream, and- in 
this case, the pulse vlll continuously grow (at least as long as 
the shock wave is in the converging part of the channel). Thus, 
smooth deceleration through the speed of sound is unstable to 
compression pulses coning from the rear. 


S'tabili'fcy of Shock Position in Channel Flows 

The previous discussion shows that, in stable channel flows, 
deceleration throu^ the speed of soimd ■will be accomplished in -a- 
shock wave.. Consider now a converging- diverging channel; There 
are, in general, two possible shock positions which ■^rt.ll yield 
equilibrium flows) that is, the shock can be either in the con- 
verging or diverging part of the channel. Experimentally, the 
shock is found to be s ■table in ■the diverging part of the channel 
and unstable in the converging part of the channel. It will be 
interesting to apply the formiilas Just developed to examine this 
problem theoretically. 

A more precise statement of the problem which will be studied 
is as follows; Consider a shock wave in a channel flow. Consider 
that a small pulse has displaced ■the shock wave by a short 
distance Since only small disturbances ■will be considered, it 

will be convenient to \ise a linearized form of the shock- velocity 
equations. If now, ■fche quantities u, a, et cetera are expanded in 
Taylor series in <!, about the equilibrium shock position, only 
the first-power terms should be ro^tained to ob^tain a linearized 
equation. Ptir'thermore , since ■the disturbances are small, the shock 
velocity can be assumed to be small compared with the local velocities 
or velocities of sound and higher powers of the shocjc velocity can 
be neglected. Performing these operations permits eqwtions (28) 
through ( 30 ) to be ■written in •the form 


(31) 


where is identical with U and ■the relaxation time t for ■the 

dt 

various reflection conditions is as follows; 
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Constant velocity: 


du, 

T = 

cbc 


2un + 3uo 
•^o o 

I - i) 


(32) 


No reflections-: 

2 S „ 2 . 
- a_ + 


<^^1.^ 3Ui^ 


6 s 


3u2o^ + ag^ (gui^ -«• 3ug^) 


dx 


1 

2 


■ ")(% * %) 


(33) 


Constant nressure: 



Values of T ^ are plotted In figure 8 against the equllihriura 

dx duj^ 

Mach number ahead of the shock. Figure B shows that *r is always 

dx 

■Dositive and the shock position is therefore stable (positive t ) 
dUi du_ 

when 2 is positive and unstable when — is negative. The 

dx dx 

shock position is thus stable in diverging channels and unstable in 
converging channels in agreement with experiment. 
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GENERAL DISCGSSICN 

Discussion of Simolifying Assumptions 


It should he pointed out attain that the test way to indicate 
the effects of the simplifying assumotions introduced in this paper 
would he to repeat sorae of the calculations made herein h'y the 
method of characteristics and to ccmuare the results obtained with 
those obtained in this patter. Since such calculations are not yet 
available, however, it would be desirable to see what information 
can be obtained frcm the results themselves in order to define their 
range of validity. This discussion is given with the central 
problem - stability of smooth transonic deceleration - in mind, but 
it will be clear that similar conclusions could be reached for cases 
where the equilibrium flow involves a shock. 

The assumptions of part I led to the trapped-pulse result and 
to the Inference that expansion waves replace the original uniformly 
decelerating flow with the steady accelerating flow appropriate to 
the channel. These assumptions have already been discussed in the 
section on -he stability of smooth deceleration through the speed 
of sound, and it has been pointed out that if this inference is 
adopted, errors due to the assumptions- introduced to facilitate the 
calculations of isentr-ouic wave propagation in part I are eliminated. 

The shock velocity formulas of appendix B are, hOTreverj only rou^ 
approximations . 

In part II, more accurate formulas for the velocity of the 
shook were found, which indicated that the trapped expansion pulses 
are consumed hy the shock and that the trapped ccmpression pulses 
continually grow as a resu-lt of the shock motion. The history of a 
ptilse has thus been divided into two Intervale. Dviring the first 
interval, the consumption of the pulse by the shock has been neglected. 
During the second interval, it has been assumed that the pulse is 
fully trapped. There will be, of course, an intermediate interval 
d^ing which the rates of the two processes are comparahle . If the 
order of magnitude of the times in which a pulse is consumed by the 
shock is much larger than the times that a ^ulse requires to approach 
this trapped state, the intermediate intei’val. will not be very important. 
The Importance of this intermediate interval can therefore be clarified 
by calculating the fraction of the pulse consumed by the shock motion 
in the time required for the trapping processes which is of the order 

of — . Calculations of this kind are presented in figure 9* 
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It can "be concluded frcsn this figure that the Intermediate interval 
does not play an .important role so long as the pulse area of the 
initial pulse is not large enough so that hi^ Mach number shocks 
are involved. 

A third group of errors was introduced in part II by applying 
the reflection conditions after the shock instead of at the end. of 
the channel. The significance of this approxlmtion can be estimated 
by comparing the results obtained in figures 7 and 8 for the various 
reflection conditicaas. The changes in reflection condition considered 
are seen to produce only a small change in the results where the 
shock Mach number is close to 1. Therefore, it seems likely that the 
results would not be greatly altered if hli^ Mach number shocks are 
not Involved by charfges in the sti-ength of disturbances in propagating 
between the shock and the end of the channel. 

The time required for disttirbances to propagate between the 
shock and the end of the channel has also been neglected and the 
significance of this anprox3.mation can be brought out by considering 
the sequence of events which would occur in the case where a very 
long channel follows abatable channel shock. Assume that a pulse 
moves upstream and disnlaces the shock. The shock then returns to 
equillbrlun at a rate which corresponds to the no-reflections case.- 
A downstream-moving mvlse will be given out by the shock while it 
is out of equilibrium, \-lhen this pulse reaches the end of the 
channe]., it will be reflected as an upscreom-movlng pulse and will 
now dlsnlac.o the shock in the onPoslte direction. (See fig. B, 
constant-pressure, or constant- velocity end conditions.) This process 
will reneat indefinitely and there is the -ooBsibility, which should 
be investigated, that a divergent oscillation could result in some 
circumstances. For the cases -vdiere this oscillation does not diverge 
(and It is clear from experiment that these cases are an Important 
group) , the application of the boundary conditions immediately 
behind the shock will suppress the damped oscillatory motion which 
would be expected for long charjaels and results in. a mcaiotonic return 
to or divergence from the equilibrium position. 


Applicatioh to Supersonic Diffusers 

An Interesting application of the present analysis will be to 
consider the maximum pulse which a supersonic diffuser flow with a 
shock in a diverging charjiel can absorb and yet retvim to its 
Initial configuration (that is, its initial shock position). Con- 
sider first the effect of expansion pulses which will move the shock 
downstream. As long as the shock is not moved beyond the end of 
the diverging part of the channel, it would be expected that the 
stability found earlier would return it to its original position. 
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The effect of ccmpressloii pulees sti'ong enougji to push the 
shock through the diffuser throat is more interesting. After the 
pulse has complete]^ interacted with the shock wave (that is, when 
the trailing edge of the pulse has just reach the shock) , the shock 
will he displaced hy an amount which can he readily estimated from 
the rule of conservation of pulse area. The results of part I 
indicate that the flow behind the shock wave will assume the alternate 
(subsonic) steady flow possible for this channel. If the pulse is 
stiff iciently strong, the displaced position of the shock will be in 
the converging part of the channel where the shock position is unstable. 
There is, of course, an unstable equilibrium position for the shock 
in the converging nart of the channel and the displaced shock will 
move away from, this unstable equilibrium position. If now the dis- 
placed shock position is d.ovnstream from the unstable equilibrium 
tiosiHon, it will move further domstreaTi and eventually return to 
the diverging part I of the charnel where it will again assume its 
stable equilibrium position. On the other hand, if the displaced 
position is upstream from the unstable equilibrium shock position, 
it will continue to move upstream and eventually convert the super- 
sonjc flow in the converging part of the channel to a subsonic flow. 

Thus, the supersonic flow with a shock in the diverging part of ilie 
channe]. is stable to compression pulses which are not sufficiently 
strong to displace the shock beyond the unstable equilibrium position 
in the converging part of the channel. 

In the practical design of supersonic diffusers, it is desirable 
to ensure that disturbances of a given magnitude (that is, a given 
pulse area) do not force the shock beyond this limiting position, 
it is, of course, also desirable that the equilibrium shock intensity 
be kept as low as possible. It appears, therefore, that diffusers 
with a long throat region, which produce a velocity distribution such 
as that shown in figure 10, should be considered. In this case, the 
compression pulse area that can be absorbed by the diffuser, that 
is, the area ABCD, is increased without increasing the shock intensity. 
However, it will be noticed that the skin friction would also be 
increased in this way. Further investigation will be necessary to 
determine the optimum throat length and shape. 


Suggestions for Future Research 

It would be interesting to perform calculations on some of the 
problems which have been studied herein by the characteristics 
method. For example, the central nroblem of this paper, stability 
of smooth deceleration throtigh the speed of sound, could probably 
be studied by this method. By starting with equations 0 ^) and ( 15 ) , 
it should be possible to trace the propagation of an upstream-movipg 
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■pulse. After a shock has heen formed in 'the pulse, it vould he 
necessary, of course, to introduce the Bahkine-Hougonolt relatione 
to trace further development. It ■would, in particular, he very 
■\raluahle to see whether more exact calculations confirm the inference 
that a tra-oped pulse converts a portion of a decelerating charnel 
flow to a steady accelerating flow. 

An experimen^bal check on part I of this theory could he obtained 
hy the use of high-speed flo^w measurements. Consider a nozzle in 
which the flow Just reaches hut does not exceed the local speed of 
sound. In such a flow, it would he expected that ■the ■trapped pulses 
(or nearly trapped pulses) would he present due to natural or 
artificial disturhances coming fi‘om ■the rear. It would he anticipated 
that ■the characteristic ■triangular shaije of ■trapped pulses (which 
would he expected also if density were plotted ins^tead of velocity) 
predicted from this ■theory would he readily dbservahle . 

An experimental Investigation of supersonic diffusers wi^th long 
throats should he conducted. This study would he particularly 
interesting in ■the case of high Mach number diffusers wi^th variable 
geometry. When the geometry is variable, the only limitation on 
the ailnimum shock in^tensity woixld he that adequate s^tahility to 
disturbances woxild he required. In many applications, a high dis- 
turbance level will he present and an effort should he made in 
experiments of this kind to simulate the actual disturhances which 
wo^uld he present in the antic i-ca ted practical application. 

As has been pointed out previously, a study of unsteady transonic 
flow in ■two or three dimensions would probably do much to clarify the 
problems of 'the formation of normal shocks in these flows. 


CCKCLtroiNG BEMAEKS 


The propagation of u-ostre am-moving pulses in the sonic region 
of a nonviscous decelerating channel flow has been studied. It is 
nointed out that a fairly complete history of a small pulse can he 
obtained by considering the two types of processes which occur 
separately (parte I and II) . 

In part I, Eiemann's theory of the propagation of finite- 
amplitude disturhances in a homogeneous medium has been extended 
to the case where upstream-moving pulses are superposed on a 
decelerating charnel flow. It is concluded ■that a pulse approaches 
a trapped state in which it cpnver^ts a portion of the channel ■fco 
the accelerating flow which is an alternate steady- s^tate solution 
for the channel. Shock waves are formed as ■the pulse is propaga^ted, 
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as in EiS'iiann's problem. If It is aesuin©d that the intei^lty of 
the shocks is very small, then, to a first approximation, it is 
shown that the shocks neither consume the pulse nor make it grow. 

It is also shown that when a pulse interacts with a weak shock 
which forms part of the original eqtiilihrium, channel flow, the 
pulse is transformed into a displacement of the shock wave . These 
facts can he summed up in the following approximate rule: the 

quantity J Qdx (called pulse area) is conserved in the propagation 

of upstream- moving pulses in the sonic region of a channel flow. 

The theory is applied to trace the development of a shock in a 
de Laval nozzle as the hack pressure is lowered. It is shown that 
a shock forms inevitably if the hack pi’essure is lowered helow the 
value which first produces sonic velocities at the throat. It is 
hoped that this calculation will he usef\aZl. in leading to the 
solution of the more difficult prohlems of the formation of shocks 
in two- oz’ tliree-dimensional-flow prohlems. 

The central problem considered in this paper is the stability 
of channel flows involving smooth (shock-free) d.eceleration through 
the speed of sound. For this problem, the analysis of part I indicate, 
tliat a pulse coiling from the rear of the channel will remain 
pen.ianently trapped In the sonic region. Thus, from part I, uhe flow 
has neutral, stability to upstream- no ving pulses. Therefore, it 
appears that a more accurate analysis must he made before stahJJ.ity 
problems can be considered. 

• In part II, the aspuription of weak shocks (used in deriving the 
rule of conservation of pu3 se area in part I) is dropped and the 
motion of shock TOves Ifi considered -more accurately. It is necessary 
to assume some type of reflection condition at the downstream channel 
end before the pi'oblem of shock motloa is competely specified. It 
is made plausible from the results that little difference in the 
shock velocity occurs when a range of reflection conditions are 
assumed. Therefore, for the exploratory study of stability problems, 
a simplified treatment of reflections la used which permits the 
calculations to be made in closed form rather iiian by the laborious 
method of characteristics. 

By combining the various reflection conditions with the 
Eankine-Hugoniot relations, the velocity of shock waves in channel 
flows which are out of equilibrium is computed. For cases where 
a trapped pulse is superposed on a smooth deceleration thi'ough the 
speed of sound, it is fomd that the shock always moves upstream. 

Thus, trapped expansion pulses are consumed by the shock motion and 
trapped compression pul.ses inevitably .grow. It is ccjncluded, there^ 
fore, that smooth ti’ansonic declaration is unstable to ccsnpreasion 
pulses coming from the rear of the channel. 
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The sta'bility to small dlstiirbances of chdimel flows Trfilch, 
at equllihrlimi, involve shock Waves is also considered in part II. 
It is shown that the shock position is mistahle in converging 
channels and is stable to small disturbances in diverging channels . 
It is made plausible that the disturbance level can affect the 
minimufa shock Intensity which can be attained in a practical super- 
sonic diffixser. The theory indicates that supersonic diffusers 
with long throats may peimit a lower shock intensity and thus have 
a hi^er efficiency when a high dist^^rbance level is present. 
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APPENDIX A 


SYMBOLS 


X distance measured along the axis of the channel; x is taken 
positive in downstream direction 


t time 



u velocity in the x direction 

a local velocity of sound 

P local' density 

Y ratio of heat capacity at constant pressure to heat capacity 
at constant volume, equal 7/5 for room-temperature air; 
this value of 7/5 is assumed throughout thj.s report 

Mq Mach number (uo/ao) 


u» r 


a* = 


U - Uq 

a - a„ 


I 


u + 5a for section on disturbances superposed on a homogeneous 
raedi-ura 


u’ + 5a’ for the section on disturbances superposed on a fldv; 
field 

[u - 5a for section on disturbances super^sed on a homogeneous 
medium 


ju* - 5a’ for the section on disturbances superposed on a flow- 
field 


5 ^ 

dt 


denotes defdvation with respect to time v/hen point under 
observation moves with local propagation velocity 
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3h 


^ displacement of a sho.ck from its equilibriimi position 
U velocity of a shock wave in stationary coordinates 
V flow velocity in coordinates moving with a shock wave 
S local channel area 



E 


du^ 

dx 


Subscripts: 

s stagnation conditions 

0 local equilibrium flow conditions 

cr equilibrium conditions for Mai 

1 conditions ahead of a shock 

2 conditions behind a shock 


ig equilibrium flow condition upstream of a 

shock at the equilibrium shock position 


2q eqviilibrium flow condition dovmstream of 
a shock at the equilibrium shock 
position 

1» equilibrium supersonic flow condition 
^ upstream of a shock calculated at 

the displaced shock position 

2,^ equilibrium subsonic flow condition 

downstream of a shock calculated at 
the displaced shock position 


Defined in 
figure 6 


A bar over a quantity denotes a lover limit of Integration 
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APPENDIX B 

APPROXIMATE FORMULAS FOR THE VELOCITY OF PROPAGATION 
OF SHOCK WAVES FOR USE IN PART I 


The expct computatioii of shock velocity, in general, involves 
cumbersoEi© algebra; therefore, for use in nart I, approximate 
expressions wil], he developed for the shock velocity. First, con- 
sider the case vhers an upstream-moving disturbance involving a 
shock is superposed on a smooth steady flow Conservation of 

energy across the shock can he written 

= Vg^ + 5fia^ (Bl) 

where the velocities are measured relative to a coordinate system 
moving with +he shock velocity U. Returning to the fonner convention 
of measuring velocities relative to stationary coordinates gives, 
for equation (Bl) 


(u;^ - U)^ + 5a^^ = (ti2 “ U)® + 582^ (B2) 

Introducing the steady flow and disturbance quantities, as before, 
gives 


Uq^ + 2u'2_Uq + Vi'j^ - 2(Uq + u'3_)U + + '^a.^ + lOa^a’^ + 58'^^ 


= Uq^ + + U^ + 5a^^ + lOa^^a'g + 

CB3) 


If only upstream-moving disturbances are present, P » 0, u' = “58', 
and equation (B3) can be simplified to 


U = u. 



8q + Qs) 


(B4) 
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In this calculation dovnstream-moving disturhances created hy 
the shock itself have heen neglected. 


The case -where a steady flow involves a shock and where sntall 
upstream-moving disturhances interact with the shock also yields a 
simple result. In this case, the steady flow velocities on the two 
Bides of the shock are different and will "be denoted hy u. 

h 

and Uo„. (See fig. 6.) Equaticn (B2) can now he written 


u, ^ + Su'nU- + - 2/ui + u'n\U + + 5eh 

X \ ^ Is 7 j 




s .?> 


- 2/uo + u'o^u + + 5ao ® + lOa^s a'c + 5a'o" 

V ^ J £ ^ 


(B5)- 


A 


Since the disturhances are now assumed small compared with the corre- 
sponding steady flow quantities, the shock velocity will also he small 
ccmpared with the steady-flow velocity. If second-order small 
quantities are neglected and it is noted that ^ ^®1 ® ^ 5^2 ^ 

equation (B5) can he written ^ S ^ ^ 




U 


- ui^u'i + ?a2^a*g - 5aj_^a»i 




- til 


(b6) 


If only upstream-moving distxurhances are assumed to he present, 



■5a and u’^ 


Equation (b 6) can then he written 
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If the variation, of the steady flow velocities with ^ is neglected, 
equation (B7) can he written 




(B8) 


Thus, in going from equation (B7) ho equation (B8) effectively, the 
effects due to stahilily of the shock wave have been neglected in 
comparison with the disturbance effects. 

■*A limiting valine of the quantities in parenthesis in equation (B8) 
as all the velocities apnroach the local speed of sound a^^ can be 

found by differentiating the steady- flow re.latlonship 

Up ® + 5a2 ^ “ 5ag^ 
o 

Thus, there is obtained, if “ ^cr ° 


®cr ^“2o * ° 



Similarly, differentiating the shock relation 



gives 


Aui = -AUo 
o o 

Prom these first approxHmatlons, it is found that, for values of M 
close to 1 
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and, similarly, 


U' 


u 


1 " 

12 22 

1 " ^ 




Thus, equation (BB) becomes 


U = ^CQi + Qg) 


(B9) 
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APPEND3X C 

CAIOTLAHONS OF THE FLOW SHOWN IN FIGURE 5 


For purposes of discussion, the. disturhances viU he divided 
into three phases: (l) a leading phase vhlch originates while the 

pressure at the end of the nozzle is fallingj (2) a steady phase 
which originates after the nozzle hack pressiire has reached a 
steady value; and (3) a shock phase which appears, in tho case 
assumed, at the Juncture of the two previous phases. 

In order to simplify the calcxilations of the leading phase, 
the rate of pressure drop (or velocity Increase) at tho nozzle end 

was chosen to make e or constant. The leading phase thus 

appears as a straight line end can he seen fresm eq.uation (2.3) to 
remain a straight line as long as the dlsturhance remains in a 
uniformly decelerating field. The motion of the leading edge of 
the dlsturhance in a short time At will he At(uQ - a^) . The end 
point of tho leading phaso can ho followed similarly when the steady 
phase propagation has been calculated. 


In the steady phase, no changes are being propagated in either 

direction; hence, a 0. Thus, a single curve will he an envelope 

for the disturhances at successive times. This curve can he found 


by integrating equation (16), while making iise of the fact that *» 
Equation (16) can then ho written 


du' 

dx 



du 

o 

dx 


u' 


1*0 ■ 


+ I u' 
5 


(Cl) 


du 

For the nozzle assumed — * is a constant which will he denoted 

dx 

by then, since M is close to 1 

D ° 


0 . 


u_ ^ a^i w ^ ? 
° ° 5 h 



4o 
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and 



g. 


This, Deviation can he written 

o£ du'*' ■+ u' d£ + hu* du' » 0 


and' Integrated to give (with quantities refeirring to the lower liait 
denoted hy "bars) 


xu'^- x5* 


2 2 


B 0 


From this equation, u' can be obtained as a function of x 
from the known valiis of. u* at the end of the nozzle. This result 
has been plotted as the envelope in figure 5* It will be noticed 
ths^t tho slope becomes infinite at the point where the steady-phase 
curve reaches the local speed of soimd. This result can be readily 
verified from equation (Cl) where the dencaninator of the right-hand 
side vanishes at the local speed of sound. This point, then, is 
where a shock first fonns. The propagation of the dlstiurbance up to 
the point where a shock forms' can now be readily calculated. Onee 
formed, the shock will connect the leading phase with the steady 
phase. 

The shock position and intensity at any time can be calculated 
step by step as follows: First, the leading phase is determined 

from the known propagation velocity of its leading edge and its 
slope as found from equation (23) . Second, the shock position is 
determined by adding to the previous shock position an increment 
corresponding to the shock velocity (equation (B4)). Third, the 
shock is then a vertical line (in fig. 5) connecting the leading 
phase with the steady phase -at the known shock position. 
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(a) Disturbance approaches shock. 
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("b) Interaction of disturbance with shock. 


Figure !•- The interaction of an expansion pulse with a nozzle 
shock. The analysis of part I shov/s that the pulse area 
which disappears on the downstream side of the shock is equal 
to the pulse area which appears due to the shock displacement 
on the upstream side of the shock. Note that the pulse areas 
plotted here are \.f Q dx since u* = S , 


Decelerating flow field 
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Accelerating flow field 






Figure 5 *“ Schematic illuatratlon of four successive positions In the 
progress of expansion and compression pulses In a uniformly decel- 
erating field. The pulse area is conserved and pulse shape approache 
an isosceles triangle. When the leading portion of the expansion 
pulse or the trailing portion of the compression pulse approaches the 
local velocity of sound, the pulse approaches a "trapped” state* 
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Exact alternate , 
equilibrium flow 


.Expansion phase 
of a trapped 
pulse E = ••l.O 


-Original equilibrium 
uniformly decelerating 
flow 
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Figure Comparison of the velocity distribution 

for the trapped expansion wave with the accelerating 
equilibrium flow for the same channel# 
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Fig. 5 




Figure 5«” development of a normal shook In a de Laval nozzle 
as t^ baok pressure is lowered at the nozzle end. The envelope 
curve is calculated in appendix C* The successive expansion wave 
positions are calculated by finding the propagation velocities of 
the end points. The time intervals between successive positions 
are equal up to tTie formation of shook and double thereafter. The 
shook velocity is calculated from the approximate equation (Bij.). 


cr 


Fig. 6 
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X 


Figure 6#“ The notation used in shook- velocity problems* 
For the case where no initial shock is Involved in 
the equilibrium flow Ut = u^ = u Explicit defi- 
ne 

nltlons of the subscripts are given in appendix A* 
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Rsflaotloa oondltlon at 
doimatream channel and 


— No reflections 
(equation (53)) 

— Constant velocity 
(equation ( 52 )) 

Constant pressure 

(equation (5^4-)) 
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Equlllbrlum-sbook Uaoh number 

Figure 8,- The relaxation time for the return of a shock to Its 
equilibrium position. The value of t is computed from the 

linearized shock motion equation ^ =0, The value 

at 

of t Is seen to be positive (stable shook position) for 
'i«lE 

a- >0, that is, for diverging channels, and negative 

dx 

(unstable shook position) for converging channels. For com- 
parison, the relaxation time for the approach of a pulse to 

the trapped stats Is ' t. liu-, shown by the dashed line. 
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Fig. 9 
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9 

Figure 9«- Fraction of the area of a trapped pulse which Is con- 
sumed In the time — • • This result Indicates that If the 

dx 


shook Intensity Is small eno\igh, only a small fraction of the 
pulse is consumed before the pulse Is fully trapped. In such 
oases, the division into trapping and consumption processes 
will be justified. 


i k 




Fig. 10 
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Figure 10,- Suggested long-throat supersonic diffuser. The 
limiting compression pulse ttiat could be absorbed by this 
diffuser would have the pulse area, ABCD, Note that the 

pulse areas plotted hero are j/'q since here u* = ^ • 
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